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Abg tract 

The  effect  of  nose  shape  on  two  flare 
stabilized  projecciles  was  studied  using  a 
Parabolized  Na v ie r -S t oke s  code.  Pressure 
coefficients,  forces  and  moments,  skin  friction 
coefficients,  and  Stanton  number  calculations  are 
presented  for  che  hemis phere-cy  l  inder-f  1  are  and 
the  cone-cy 1 inder-f 1  re  configurations.  Pitching 
moment  and  static  margin  plocs  versus  flare  angles 
and  flare  lengths  are  presented  in  a  parametric 
study  to  show  aerodynamic  stability  effects. 
Distinctly  different  flow  field  values  of 
pressure,  local  Mach  number,  and  dynamic  pressure 
were  generated  by  the  two  different  nose 
configurations.  These  flow  field  values  just 
upstream  of  the  flare  will  be  examined.  The  code 
demonstrated  its  value  as  a  design  tool  by  making 
a  clear  distinction  between  aerodynamica l ly  stable 
characteristics  for  this  variety  of  nose  shapes, 
flare  angles,  and  flare  lengths. 

Introduction 

The  stabilizing  effect  of  flared  afterbodies 
has  long  been  incorporated  in  the  design  of 
rockets  and  projectiles.  For  short  projectiles  at 
relatively  high  Mach  numbers,  the  flare  angle  and 
the  flare  length  can  play  a  significant  role  in 
the  projectile's  aerodynamic  stability.  Extensive 
experimental  studies  have  been  made  on  cone- 
cylinder-flared  projectiles  at  Mach  numbers 
ranging  from  2  to  41_J  and  a  Parabolized  Navier- 
Stokes  (PNS)  code  has  been  used  for  numerical 
■studies.*1’3  Most  of  these  studies  involved 
conical  nosed  bodies  at  lower  supersonic  speeds. 

This  paper  will  present  aerodynamic 
characteristics  determined  numerically  for  two 
flare  stabilized  configurations  at  Mach  7  and  at 
2*  angle  of  attack;  a  hemisphere-cylinder-flare 
and  a  cone-cylinder-flare.  Mach  7  was  chosen 
because  it  was  high  enough  to  yield  characteris- 
ticly  high  Mach  number  effects  yet  low  enough 
that  real  gas  effects  can  be  ignored.  The  low 
angle  of  attack  was  selected  to  generate  pitching 
moment  data.  The  effect  of  the  nose  shape  on  each 
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projectile's  aerodynamic  stability  will  be  consi¬ 
dered.  Within  that  context,  the  flare  angle  and 
the  flare  length  will  be  varied  to  study  the  effect 
of  these  two  parameters.  A  brief  description  of 
the  PNS  code,  code  inputs,  configurations,  and 
flow  conditions  will  be  presented  in  the  second 
section.  Aerodynamic  results  will  be  compared  for 
each  configuration  in  the  third  section.  Finally, 
conclusions  to  be  drawn  will  be  presented  in  the 
last  section. 

PHS  Code  and  Configurations 

A  Parabolized  Navier-Stokes  code^  used  for 
this  numerical  study  generates  finite  difference 
approximations  to  steady,  three-dimensional 
solutions  of  supersonic  flow  over  arbitrarily 
shaped  bodies  at  high  Reynolds  number  as  long  as 
there  are  no  large  subsonic  or  axially  separated 
regions.  From  an  initial  condition  starting  solu¬ 
tion,  the  code  marches  downstream  from  any  given 
axial  location  on  the  body.  For  the  blunt  hemis¬ 
pherical  nose  configuration  in  Figure  1(a),  the 
code  needed  an  input  starting  solution  obtained 
from  a  thin-layer  unsteady  Navier-Stokes  code- 
For  the  sharp  conical  nose  configuration  in  Figure 
1(b),  the  PNS  code  has  a  self-start  capability 
wherein  initial  starting  planes  are  calculated 
within  the  code. 

The  PNS  code  is  a  space  marching  code  which 
uses  a  numerical  method  originally  developed  at 
NASA  Ames  Research  Center  by  Schiff  and  Steger. 
This  technique  obtains  finite  difference  solutions 
by  solving  implicit  difference  equations 
noniteratively  by  way  of  local  linearization  of 
the  flux  vector.  The  resulting  system  of 
algebraic  equations  are  solved  using  the  Beam- 
Warming  method”  of  implementing  the  alternating 
direction  implicit  scheme  using  the  delta 
formulation.  Several  additional  modifications  to 
the  code  have  been  made.^®  * 

The  PNS  flow  field  predictions  for  the 
hemisphere -cy  1  i  nde  r  -  f  1  a  re  were  obtained  by 
marching  from  a  starting  solution  generated  by  a 
thin-layer  unsteady  Navier-Stokes  code.  This 
starting  solution  was  obtained  at  an  axial 
location  at  che  hemisphere-cylinder  shoulder  with 
an  algebraic  grid  of  19  by  30  points  in  the 
circumferential  and  radial  directions, 
respectively.  This  data  was  then  interpolated  to 
obtain  an  algebraic  grid  of  19  by  45  points.  The 
data  was  obtained  at  0*  angle  of  attack  and  then 
rotated  co  a  2*  angle  of  attack.  For  a 
hemispherical  nose,  this  rotation  presents  no 
problems  theoretically.  Thirty-five  of  the  45 
radial  grid  points  were  clustered  near  the  wall  to 
resolve  the  boundary  layer.  A  marching  step  size 
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of  0.01  cylinder  diameters  was  used  up' to  a  loca¬ 
tion  one-half  cylinder  diameter  aft  of  the 
hemisphere-cylinder  shoulder.  This  small  step 
size  captured  the  strong  gradients  in  Che  flow 
field  in  this  shoulder  region.  The  step  size  was 
then  increased  to  0.03  diameters  for  marching 
along  the  cylinder  co  a  plane  0.13  diameters 
upstream  of  the  cylinder-flare  juncture.  A 
reduced  step  size  of  0.02  diamecers  was  used  to 
adjust  for  che  rapid  pressure  rise  in  the  region 
at  che  flare.  This  step  size  was  close  to  one- 
third  the  boundary  layer  thickness  predicted  by 
che  PNS  code  at  this  point.  Analysis  by  Viegas 
and  Horstman  indicate  that  a  step  size  one-half 
the  boundary  layer  thickness  is  sufficient  to 
resolve  two-dimensional  shock/boundary  layer 
interactions.14  Calculations  were  made  on  0,  5, 
10,  15,  20,  and  25  degree  flare  angles.  Attempts 
to  calculate  a  30”  flare  angle  failed  with  axial 
flow  separation  characteristics.  This  is 
consistent  with  the  PNS  code's  ability  co  detect 
the  onset  of  incipient  separation  as  evaluated  by 
Coctrell  and  Chapman."*  This  failure  may  also  be 
attributed  to  the  3eam-Warming  type  of  algorithm 
which  tends  co  introduce  errors  in  the  form  of 
local  flow  property  oscillations  as  a  result  of 
the  cencral-differencing  of  fluxes  across  discon¬ 
tinuities.1-*  Versions  of  the  PNS  code  using  an 
upwind-differencing  scheme  to  alleviate  this 
problem  are  being  employed.  The  code's  smoothing 
terms  were  as  follows:  explicit  damping  parameter 
of  0.2,  implicit  damping  parameter  of  0.4,  sta¬ 
bility  coefficients  of  0.25.  The  nominal  v^ilue  of 
the  grid  spacing  at  the  wall  was  0.5  x  lO"3. 


The  PNS  flow  field  prediction  for  Che  cone- 
cylinder-flare  configuration  was  obtained  by 
marching  from  a  solution  generated  by  the  code's 
self-start  option  for  sharp  conical  noses.  This 
starting  solution  was  obtained  at  an  axial 
location  1  cylinder  diameter  from  the  nosetip  at  a 
2”  angle  of  attack.  The  algebraic  grid  was 
19  by  45  points  in  the  circumferential  and  radial 
directions,  respectively.  Thirty-five  of  the  45 
radial  grid  points  were  clustered  near  Che  wall  to 
resolve  the  boundary  layer.  This  grid  is  set  up 
exactly  the  same  as  Che  grid  for  the  hemisphere- 
cylinder-flare.  However,  the  volume  of  Che  grid 
for  the  cone-cylinder-flare  case  is  much  smaller 
than  the  volume  of  the  grid  for  Che  hemisphere- 
cy 1 inder-f lare  case  due  to  che  smaller 
distance  of  che  bow  shock.  Hence,  the  grid  points 
clustered  near  the  wall  for  the  cone-cv 1 inder- 
flare  are  packed  closer  together.  A  marching  seep 
size  of  0.02  diameters  was  used  co  a  plane  0.8 
cylinder  diameters  aft  of  the  cone-cylinder 
shoulder.  This  small  step  size  resolved  the  large 
gradients  in  the  flow  field  as  the  solution 
marched  up  over  the  shoulder.  The  step  size  was 
then  increased  to  0.04  diameters  for  marching 
along  the  cylinder  to  a  plane  0.14  diameters 
upstream  of  the  cylinder-flare  juncture.  The  step 
size  was  reduced  back  down  to  0.02  diameters  to 
adequately  resolve  che  rapid  pressure  rise  in  Che 
flare  region.  For  Che  20*  flare,  Che  step  size 
was  kept  at  0.04  diameters.  A  flare  angle  of  25* 
was  attempted  but  che  code  indicated  che  onset  of 
axial  flow  separation  just  as  Che  30*  flare  angle 
attempt  for  the  hemisphere-cylinder-flare  did  and 
the  reasons  cited  at  chat  time  are  applicable 
here,  also.  The  smoothing  ceres  were  kept  the 
same  as  those  used  in  the  hemisphere-cylinder- 
flare  study  for  all  cases. 


Hemisphere -Cylinder-Flare 


Figure  2  shows  che  shock  structure  for  the 
hemi sphere-cy l inde r- 1 5  degree  flare  sc  Mach  7. 
The  outer  boundary  shock  was  fit  by  the  unsteady 
Navier-Scokes  code  for  the  hemisphere  nose  segment 
and  che  PNS  code  fit  the  shock  for  che  cylinder- 
flare  segment.  Note  the  small  discontinuity  where 
Che  shock  fits  are  joined  together  at  tr.e  shoulder 
region.  This  is  not  unusual  when  shifting  from  an 
unsteady  Na vier-Stokes  code  to  a  PNS  code.  The 
unsteady  code  captured  the  subsonic  region  of  the 
nose.  The  sonic  lines  were  estimated  from 
locations  where  the  total  Mach  equaled  1.  The 
embedded  secondary  shock  locations  were  captured 
by  the  PNS  code  and  the  estimated  locations  in  the 
flare  region  were  obtained  from  P  /  P  £  n  f 
distributions  beeween  che  body  and  the  shock  at 
the  leeward  and  at  the  windward  planes.  These 
captured  shocks  are  smeared  across  grid  points  and 
sharp  discontinuities  are  not  apparent.  Figure  3 
shows  pressure  coefficient  distributions  at  the 
surface  at  three  circumf erent ia 1  stations:  a 

windward  ray,  a  90*  or  side  ray,  and  a  leeward 
ray.  As  anticipated  at  this  small  angle  of 
attack,  there  is  a  slightly  larger  Z0  on  the 
windward  side  which  decreases  slowly  as  one  moves 
circumferentially  to  the  lower  on  tr.e  leeward 
side.  From  a  C  of  1.83  at  the  nose  stagnation 
poinc,  the  flowHfield  pressure  undergoes  a  rapid 
expansion  over  the  nose  shoulder,  decreases 
gradually  along  the  cylinder  body,  and  jumps 
quickly  at  the  cylinder-flare  juncture.  Note  once 
again  that  the  joining  of  unsteady  Navier-Stokes 
code  data  with  PNS  code  daca  results  in  a  slight 
C  discontinuity  in  the  shoulder  region.  Figures 
4-6  show  pressure  coefficient  distributions  at  the 
same  three  circumferential  stations  on  the 
cylinder-flare  segment  for  flare  angles  ac  0  and 
15,  5  and  20,  and  10  and  25  degrees.  Note  how  the 
pressure  jumps  at  the  cylinder-flare  juncture. 
For  lower  flare  angles,  the  pressure  jumps  and 
then  decays  as  expected  for  uniform  flow 
approaching  the  flare  (i.e.,  jump  to  pressure 
associated  with  a  2-D  wedge  and  then  drop  down  to 
pressure  associated  with  a  cone).  At  the  20  and 
25  degree  flare  angles,  there  appears  a  deviation 
from  this  trend.  This  is  due  to  nonuniform  flow 
approaching  Che  flare  and  will  be  discussed 
further  in  the  comparison  with  the  cone-cylinder 
case.  Attempts  to  calculate  a  30*  flare  case 
failed  with  axial  flow  separation  characteristics 
for  reasons  discussed  earlier. 

It  should  be  noted  that  this  marching  code 
cannot  predict  the  pressures  on  the  base  of  the 
flare.  If  the  pressure  distribution  on  the  base 
is  uniform,  then  the  base  pressure  effect  on 
pitching  moment  is  negligible  for  chis  small  angle 
of  attack.  Pitching  moment  curves  versus  flare 
lengths  for  various  flare  angles  are  shown  in 
Figure  7.  Stability  grows  with  increases  in  both 
flare  length  and  flare  angle.  These  results  are 
similar  to  the  stability  trends  reported  by 
Cottrell  and  Chapman. 

Static  margin  curves  were  constructed  by 
subtracting  the  center  of  gravity  (eg)  location 
from  the  center  of  pressure  locacion.  A 
homogeneous  body  was  assumed  for  the  eg 
calculation.  The  center  of  pressure  locations 


were  calculated  by  the  PNS  code.  Both  locations 
were  measured  positively  from  the  nose.  When  the 
center  of  pressure  is  farther  away  from  the  nose 
than  the  eg,  the  nose  will  pitch  downward  which 
indicates  aerodynamic  stability.  In  Figure  8,  the 
static  margin  for  0,  5,  10,  15,  20,  and  25  degree 
flare  angles  is  plotted  for  various  flare  lengths. 
For  the  0,  5,  and  10  degree  flare  cases,  these 
curves  indicate  unstable  aerodynamics  for  whatever 
flare  length  chosen.  The  10*  flare  indicates  an 
impasse  where  an  increase  in  flare  lengch  moves 
the  eg  and  the  center  of  pressure  the  same 
distance  with  no  further  stability  deterioration 
or  increase.  Noce  that  the  15,  20,  and  25  degree 
flare  cases  go  stable  rather  quickly  with 
increasing  flare  angle  and  flare  length  in  a 
manner  similar  to  that  reported  by  Cottrell  and 
Chapman.4  These  last  three  flare  angles  all  go 
stable  for  flare  lengths  over  one-third  cylinder 
diameter  long. 


Figures  9  and  10  show  typical  skin  friction 
coefficient  and  heat  transfer  (Stanton  number) 
distributions.  In  Figure  9,  skin  friction  decays 
smoothly  with  axial  position  except  near  the 
cylinder-flare  juncture.  The  expansion  of  the 
flow  at  the  cone-cylinder  shoulder  increases  the 
boundary  layer  thickness  resulting  in  a  rapid 
decrease  in  skin  friction  coefficient  which 
sceadies  to  a  smooth  decay  down  the  cylinder 
segmenc.  The  spike  in  the  skin  friction 
coefficient  at  the  cylinder-flare  juncture  is 
physically  unreal.  It  may  be  attributed  to  the 
nature  of  the  PNS  type  of  algorithm. 
In  Figure  10,  Stanton  number,  or  heat  transfer, 
also  decreases  smoothly  with  axial  position  except 
near  the  cylinder-flare  junccure.  The  spike  at 
that  junccure  is  physically  unreal. 


Cone -Cvlinder -Flare 


Figure  11  shows  the  shock  structure  for  the 
cone-cy 1 inder- 1 5  degree  flare  at  Mach  7.  The 
attached  shock  at  the  nose  and  the  outer  shock 
location  were  estimated  by  the  PNS  code  in  the 
step-back  mode  up  to  a  location  one  cylinder 
diameter  from  the  nosetip.  Planes  of  data  were 
stored  there.  The  PNS  code  picked  this  data  up 
and  marched  down  the  rest  of  the  configuration. 
Once  again,  the  estimated  embedded  shock  locations 
on  the  flare  were  approximated  from  the  P/P^n£ 
distributions  between  the  body  and  shock  at  the 
windward  and  leeward  planes.  Figure  12  shows 
pressure  coefficient  distributions  along  the 
configuration  at  the  surface  for  the  15*  flare 
case.  Note  the  Prandtl  -  Meyer  expansion  at  the 
cone-cylinder  shoulder  which  drops  pressure 
rapidly  from  a  stagnation  C  of  2.13. 
Then  there  is  a  gradual  expansion  back  to  near 
freestream  pressure  values  down  the  cylinder 
segmenc.  Figures  13-15  show  pressure  coefficient 
distributions  on  the  cy l inder-f  1  are  segment  for 
flare  angles  of  0  and  15,  5  and  20,  and  10 
degrees.  Note  how  the  pressure  data  jumps  at  the 
cylinder-flare  junccure  for  all  cases  followed  by 
a  continual  pressure  rise  normally  associated  with 
nonuniform  flow.  This  nonuniform  aspect  of  the 
flow  will  be  diacuased  in  the  comparison  section 
chat  follows.  Attempts  to  calculate  a  25*  flare 
failed  with  axial  flow  separation  characteristics 
consistent  with  the  results  found  by  Cottrell  and 
Chapman4  but  may  also  be  attributed  to  Che  Beam- 
Warming  algorithm  mentioned  previously.  Pitching 


moment  curves  versus  flare  lengths  for  various 
flare  angles  are  shown  in  Figure  16.  As  expected, 
Che  pitching  moment  grows  as  flare  angle  and  flare 
lengch  increases. 

Static  margin  curves  were  constructed  from  the 
PNS  code  data  and  are  shown  for  flare  angles  0,  5, 
10,  15,  and  20  degrees  in  Figure  17.  Note  it 
appears  chat  the  0  and  5  degree  flare  angles  will 
never  yield  aerodynamic  stability  by  increasing 
the  flare  length.  Note,  too,  that  the  10*  flare 
goes  stable  at  a  flare  lengch  approximately  one- 
third  cylinder  diameter.  The  15  and  20  degree 
flares  yield  stability  rather  quickly  as  the  flare 
lengch  increases. 

Figures  18  and  19  show  skin  friccion 
coefficient  and  heat  transfer  (Stanton  number) 
distributions.  The  spike  at  the  flare  for  both  of 
these  plots  is  not  realistic  and  may  be  attri¬ 
buted  to  flow  property  oscillation  problems 
mentioned  earlier. 


Comparison  of  Configuration  Results 

A  comparison  of  the  aerodynamic  results  must 
start  with  an  examination  of  the  outer  shock 
boundaries  in  Figures  2  and  11.  The  hemisphere 
nose  generates  a  decached  bow  shock  with  a 
characteristic  subsonic  region  while  the  cone  nose 
has  an  attached  shock.  The  difference  between 
these  bow  shock  waves  lead  to  very  different  flow 
conditions  in  front  of  the  flares.  Figures  20, 
21,  and  22  show  the  radial  distribution  of  ?/P£nf> 
local  Mach  number,  and  dynamic  pressure  (q /q f 3, 
respectively.  These  distributions  are  taken  at  a 
location  just  upstream  of  the  cylinder-flare 
juncture  for  the  hemisphere-cylinder-flare  and 
for  the  cone-cylinder-flare  configurations.  Also 
shown  in  these  figures  are  the  effective  height  of 
a  10*  and  a  20*  flare  at  the  base.  This  gives  a 
rough  idea  of  the  flow  encountered  by  these 
flares.  Figures  23(a)  and  23(b)  show  the 
algebraic  grid  at  these  locations  for  the 
he  mi  s  phe  r  e  -c  y  l  i  nde  r  -  f  1  ar  e  and  for  the  cone- 
cylinder-flare,  respectively. 


The  P/pinf  plot  in  Figure  20  reveals  the  large 
region  of  constant  pressure  each  configuration 
has.  Pressures  just  before  the  flare  are  close  to 
50Z  higher  for  the  hemispherical  nose  configura¬ 
tion.  The  strong  bow  shocks  on  each  of  the  two 
lead  to  significant  entropy  losses.  The  blunc 
hemisphere  body  has  the  greater  losses.  This  is 
illustrated  by  the  lower  local  Mach  number  and  the 
lower  dynamic  pressure  in  front  of  the  flare  for 
this  blunt  nose  compared  to  the  cone  nose  as  shown 
in  Figures  21  and  22.  The  significant  reduction 
in  dynami'  pressure  accounts  for  the  much  lower 
scatic  margins  of  the  hemi  s  phere-cy  1  inder-f  lare 
compared  to  the  cone-cylinder-flare  configuration. 


Although  the  hemispherical  nosed  body  has  the 
greater  entropy  losses  of  Che  two,  the  conical 
nosed  body  has  more  of  a  variation  in  local  Mach 
number  radially  over  the  distance  comparable  to 
the  flare  height.  It  is  this  variation  chat 
accounts  for  the  rather  smooch  pressure  coeffi¬ 
cient  rise  along  the  flare  on  the  con ica l -nosed 
body  seen  in  Figures  12-15.  This  variation  in 
local  Mach  number  is  only  present  in  a  small 
radial  region  for  the  hemispherica l -nose  body. 
Hence,  only  the  smaller  flare  angles  for  the 
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heoispherica 1 -nose  body  have  smooth  pressure  rises 
on  the  flare  as  seen  in  Figures  3-6. 


8.  Schiff,  L.B.  and  Steger,  J.L.,  "Numerical 
Simulation  of  Steady  Supersonic  Viscous  Flow," 
AIAA  Paper  79-0130,  January  1979. 


Remarks 


Based  on  the  predictions  made  with  the  PNS 
code  on  the  hemisphere-cylinder-flare  at  various 
flare  angles  and  flare  lengchs  with  flow 
conditions  of  Mach  7  and  at  2*  angle  of  attack, 
toe  following  conclusions  can  be  drawn. 


At  this  hypersonic  speed,  the  nose  shape  has  a 
distinct  effect  on  the  aerodynamic  stabilizing 
effect  of  flared  afterbodies.  Of  the  two  configu¬ 
rations  examined,  the  conica 1 -nosed  body  has 
greater  stability  and,  therefore,  more  utility  at 
hypersonic  applications.  It  appears  an  attached 
shock  will  keep  encropy  losses  to  a  minimum.  This 
provides  the  stabilizing  control  surfaces  within 
the  outer  shock  enveloping  the  body  a  more 
effective  role  in  the  projectile's  hypersonic 
mission.  Furthermore,  the  increased  stability 
obtained  at  larger  flare  angles  and  greater  flare 
lengths  for  both  configurations  have  been 
presented. 
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r  =  0.5d  /cn  m  1 .2d 

4>{  =  0°  to  25°  /cy  -  2.13d 

<t,c  «  22.50°  /f  -  0.19,  0.51,  0.83d 

Figure  1.  (a)  Hemisphere  -  Cylinder  -  Flare  Geometry 

(b)  Cone  -  Cylinder  -  Flare  Geometry 


